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ABSTRACT 
> 

Cn , Aims. We investigate the role of the luminous infrared galaxy (LIRG) and ultra-luminous infrared galaxy (ULIRG) phases in the evolution of 
^^ ' A'j-selected galaxies and, in particular. Extremely Red Galaxies (ERGs). 

Methods. With this aim, we compare the properties of a sample of 2905 A", < 21.5 (Vega mag) galaxies in the GOODS/CDFS with the sub- 
sample of those 696 sources which are detected at 24 /jm . 

Results. We find that LIRGs constitute 30% of the galaxies with stellar mass M > 1 x 10" Mq assembled at redshift z = 0.5. A minimum 
Q I of 65% of the galaxies with M > 2.5 x 10" Mq at z « 2 - 3 are ULIRGs at those redshifts. 60% of the ULIRGs in our sample have the 
characteristic colours of ERGs. Conversely, 40% of the ERGs with stellar mass M > 1.3 x 10" M© at 1.5 < z < 2.0 and a minimum of 52% 
of those with the same mass cut at 2.0 < z < 3.0 are ULIRGs. The average optical/near-IR properties of the massive ERGs at similar redshifts 
that are identified with ULIRGs and that are not have basically no difference, suggesting that both populations contain the same kind of objects 
Q ' in difi^erent phases of their lives. 
}—{ ' Conclusions. LIRGs and ULIRGs have an important role in galaxy evolution and mass assembly, and, although they are only able to trace a 
C/J fraction of the massive (M > 1 x 10" Mq) galaxies present in the Universe at a given time, this fraction becomes very significant ( > 50%) at 
^ redshifts z > 2. 
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1 . Introduction mid and far-IR backgrounds (Papovich et al. 120041 Le Floc'h 

et al. 2005"; Perez-Gonzalez et al. "^OOS" Caputi et al. 2006b 

The study of the evolution of near-infrared (IR)-selected galax- ^^. ^^^^ ^^ ^j _2006^ ^j^^ ^^^^^^.^ ^^ ^^^ ^^^j^^.^^ ^^ ^^_ 

ies is probably the most efficient way to trace the build up of yr, i ■ vu j ut^ j * j » i • * 

^ -^ ^ ^ IR galaxies with redshirt provides fundamental information to 

stellar mass with redshift (e.g. Dickinson et al. .2003. Fontana j » ju . * .■ j »• * j j 

, , ^ =■, , ' ' understand how star-formation and quasar activity proceeded 

etaL|2004|Glazebrooketal.|2004l Caputi etaLEOOa 2006a through cosmic time. 
JC06aH . In contrast to other wavelength surveys, near-IR ob- 
servations are able to detect systems with different dust con- It is so far not clear, however, what the role of bright mid- 
tent and star-formation histories. Thus, they appear as a viable IR-selected galaxies is within the history of steUar mass as- 
method to produce a relatively unbiased census of galaxy pop- sembly. The aim of this work is to assess the significance of 
ulations, from which the history of stellar mass assembly can *is role, through the study of mid-IR-selected galaxies in the 
be investigated context of the evolution of the Ks-hand galaxy population. With 
The Spitzer Space Telescope (Werner et al. 2004i is making *is purpose, we focus on the analysis of two main topics: 1) the 
possible the mapping of the mid-IR Universe from low to high "'^"iber densities of massive luminous mid-IR-selected galax- 
redshifts. Recent studies have exploited this capability to set i^^, in comparison to the total number densities of massive 
tight constraints on the intensity and the composition of the g^l^^i^s in the K,-hmd population, at different redshifts, and 

2) the relation between the luminous IR and the Extremely Red 

Send offprint requests to: K.I. Caputi; e-mail: kcaputi@ias.u-psud.fr Galaxy (ERG) phases. We adopt throughout a cosmology with 



Caput! et al.: The role of LIRGs and ULIRGs 






0.5 



Fig. 1. The percentage of Ks < 21.5 mag galaxies which are 
detected in the deep MIPS/GTO 24 //m catalogue, as a func- 
tion of redshift. The shaded and empty histograms indicate the 
percentages of galaxies with 5v(24//m) > 83 //Jy and the total 
of 24 yum -detected galaxies, respectively. The faintest galaxies 
considered here have 5'v(24yum) ~ 20;uJy, but galaxies with 
5v(24//m) < 60 (83)juJy only constitute ~ 1 1 (30)% of the to- 
tal sample. 
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= 0.3 and Qa = 0.7. A Salpeter 
(0.1 - 100) Mo is assumed. 



IC06al selected a sample of 2905 Ks < 21.5 (Vega mag) galax- 
ies in 131 arcmin^ of the Great Observatories Origins Deep 
Survey (GOODS) Chandra Deep Field South (CDFS). This 
area benefits from deep J and /T^-band coverage by the Infrared 
Spectrometer and Array Camera (ISAAC) on the Very Large 
Telescope (VLT). Other deep multiwavelength data also ex- 
ist in this field; in particular, B, V, i and z-band imaging from 
the Advanced Camera for Surveys (ACS) on-board the Hubble 
Space Telescope (HST), and 3.6 to 8.0 /zm imaging from the 
Infrared Array Camera (IRAC) on board Spitzer. 

IC06al constructed an optimized redshift catalogue for 
the total Ks < 21.5 mag sample, incorporating the avail- 
able spectroscopic (e.g. Le Fevre et al. 2004; Vanzella et 
al. l20U5t and COMBO 17 (Wolf et al. IMMl redshifts for the 
CDFS. The remaining redshifts were obtained with the public 
codes HYPERZ (Bolzonella, Miralles & Pello 20001 and BPZ 
(Benitez l2000> . using multiwavelength photometry. The com- 
parison of the redshifts for the sources with spectroscopic data 
showed that the photometric estimates had very good quality, 
with a median of \zspec - Z;Vm/I/(1 + Zspec) - 0.02. The code 
HYPERZ has simultaneously been used to model the spectral 
energy distribution (SED) of each source, with the Calzetti et 
al. (I2000> reddening law to take into account dust-corrections. 



The Ks < 21.5 mag sample is cleaned of galactic stars. 
It also excludes those X-ray sources whose SEDs could not 
be satisfactorily fitted with the HYPERZ stellar templates, as 
their optical/near-IR light could be contaminated by an active 
(AGN/QSO) component. For these sources, no reliable mass 
estimate can be obtained from the modelled SED. The rejected 
X-ray sources constitute ~ 2% of the total Ks < 21.5 mag 
sample. We refer the reader to lC06al for further details on the 
selection and analysis of the Ks < 21.5 mag galaxy sample. 

Observations of the CDFS have also been carried out with 
the Multiband Imaging Photometer for Spitzer (MIPS; Rieke 
et al. .2004) . as part of the Guaranteed Time Observers (GTO) 
program. The MIPS GTO lAjim catalogue achieves 80% com- 
pleteness at a flux limit of 5v(24yum) ^ 83/vJy (see details in 
Pap oviche tal. .20043 . 

IC06bl cross-correlated the IC06al sample (and all the ex- 
cluded Ks < 21.5 mag stars and AGN/QSO) with the deep 
24 jum catalogue for the CDFS. ~ 94% of the 24 jum objects 
with flux 5'v(24//m) > 83 //Jy were identified with a 
Ks < 21.5 mag counterpart. This allowed for the charac- 
terization of nearly all the brightest sources composing the 
24 jum background in the GOODS/CDFS. 

The fraction of 24 yum galaxies with multiple Ks-h?inA asso- 
ciations is small. C06b found that, within a matching radius of 
2 arcsec, < 8% (0.3%) of the 24 jum galaxies could be associ- 
ated with two (three) different ^5-band galaxies. On the other 
hand, ~ 95% of the associations can be done restricting the 
matching radius to 1.5 arcsec. Using this smaller matching ra- 
dius, the number of double identifications is only < 3%. In all 
cases, C06b.considered that the closest /T^-band source to each 
24 jum galaxy was the real counterpart. 

In addition, using empirical relations between the mid-IR 
and bolometric IR luminosities (Chary & Elbaz 2001 Elbaz 
et al. I2002> . IC06bl obtained IR luminosity (L/r) estimates and 
derived star-formation rates (SFR) for most of their galaxies. 
In agreement with other works (Le Floc'h et al. 2005 Perez- 
Gonzalez et al. 12005 > . they found that the mid-IR output was 
dominated by luminous IR galaxies (LIRGs, 10'' L© < L < 
10'^ Lq) at redshift z ~ 1 and by ultra-luminous IR galaxies 
(ULIRGs, L > 10'2 Lq) at redshift z > 1.5 - 2.0. 

In this work we analyze the role of the sub-sample of the 
Ks < 21.5 mag galaxies identified with a 24 jum counterpart, 
within the total Ks < 21.5 mag galaxy population. We quantify 
the importance of galaxies with different IR luminosities in the 
evolution of massive systems at different redshifts. The sub- 
sample of 24 /vm -detected galaxies contains 696 out of 2905 
Ks < 21.5 mag galaxies. 

The percentage of the Ks < 21.5 mag galaxies which 
are 24 //m -detected as a function of redshift is shown in 
Figure n The shaded and empty histograms indicate the per- 
centages of galaxies with 5'v(24jum) > 83juJy and the to- 
tal of the 24 yum -detected galaxies, respectively. The faintest 
sources included in our catalogue have fluxes SyilApm) ~ 
20juJy, but galaxies with 5v(24/vm) < 60(83)yuJy only con- 
stitute ~ 1 1 (30)% of the total sample. The number of spu- 
rious sources in the 24 //m catalogue becomes non-negligible 
at fluxes 5v(24jum) < 83 //Jy and starts to be important at 
Sv{24fj.m) < 60juJy (Papovich et al. 2004). Nevertheless, we 



Caputi et al.: The role of LIRGs and ULIRGs 



note that the reliability of all the 24 //m sources analyzed here is 
assured by our Kg and other band detections. Even so, there still 
could be the possibility of a false 24 jum source which is asso- 
ciated just by chance with a K^ counterpart. However, we men- 
tioned above that the fraction of 24 //m sources with a double 
association within 1.5 arcsec is < 3%. We can consider this per- 
centage to be representative of the maximum fraction of false 
identifications. Thus, we estimate that the fraction of spurious 
Sv{24fim) < 60yuJy sources associated just by chance with a 
Ks galaxy must constitute < 0.03 x 11%=0.33% of our total 
sample. Then, those plausible spurious sources do not consti- 
tute a concern in this work. 

We see in Figure [l] that ~ 20 - 25% of our K, < 21.5 
mag galaxies at redshifts z < 1.5 are 24 jum -detected at the 
depth of the MIPS/GTO CDFS catalogue. The percentage of 
24 /im -detected galaxies increases at higher redshifts, with a 
maximum of ~ 43% at z ~ 2.0 - 2.5. The magnitude limit of 
our Ks-hand survey imposes that we only see massive galax- 
ies at high redshifts. Above redshifts z = 2 and z - 3, we 
are only strictly complete for stellar masses M > 6 x 10''^ M© 
and M > 1.3 X 10" Mq, respectively (cf. C06a and Figurel^. 
Thus, the increasing percentage of IR galaxies with redshift in 
particular implies that the fraction of massive galaxies with IR- 
activity has been significantly larger in the past. We further dis- 
cuss this issue in Section |3 However, it should also be noted 
that this positive selection effect on high-redshift galaxies is 
very probably favoured by the presence of policyclic aromatic 
hydrocarbon (R\H) emission features entering the 24 /im -filter 
passband (cf. "C06b"i. PAH emission lines characterize the in- 
terstellar medium of star-forming regions (Desert, Boulanger 



& Puget [T990b . The presence of PAHs implies that a substan- 
tial fraction of the Ks < 21.5 mag galaxies at z ~ 2.0-2.5 
have IR activity mainly due to star-formation. PAH emission in 
high-z galaxies has been spectroscopically confirmed in dififer- 
ent recent works (e.g. Houck et al. 120051 Yan et al. 120051 Lutz 
et al. ;2005) . 

3. The normal- versus active-galaxy separation 

As we mentioned above, those X-ray sources with no satis- 
factory modelled SED have been left out of the Ks < 21.5 
mag sample. This implies that 40 X-ray-detected 24 yum sources 
(including X-ray AGN/QSO and other galaxies with undeter- 
mined X-ray classification) are excluded from the sub-sample 
analyzed here. However, we note that our 24 /vm sub-sample 
still contains 26 sources X-ray classified as AGN/QSO (with 
satisfactory HYPERZ SED fits). 

To test the completeness of the X-ray detections to sep- 
arate all the 24 //m -detected active galaxies, we construct a 
colour-colour diagram based on different IRAC passbands for 
the high redshift (z > 1.5) sources. The separation of nor- 
mal and active galaxies in such a diagram should be produced 
by the mid-IR excess characterising the latter type of objects. 
At lower redshifts, an AGN/QSO separation based only on 
the X-ray classification is expected to be more complete than 
at higher redshifts. On the other hand, the existence of star- 
forming galaxies with PAH emission lines within the IRAC 
passbands complicates the colour-colour separation at low red- 
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Fig. 2. IRAC-based colour-colour diagram for the Ks < 21.5 
mag sources with redshifts z > 1.5. In this plot, we label as 
'normal' to all those galaxies which are not X-ray classified as 
AGN/QSO. Galactic stars have also been added for a compari- 
son. 

shifts, as these emission lines might mimic the AGN IR colour 
excesses. Thus, we restrict the analysis of the colour-colour di- 
agram to sources with z > 1.5. Figure |5] shows the location 
in this diagram of different types of 24 /vm sources and all the 
Ks < 2 1 .5 mag sources, at z > 1 .5 . We clearly see from this di- 
agram that the majority (~ 70%) of the 24 /im -detected X-ray 
AGN/QSO lie in a segregated region, typically with colours 
[5.8] - [8.0/im] > 0.2 (AB mag). A few apparently normal (i.e. 
non X-ray classified AGN/QSO) galaxies appear within this 
region. The vast majority of normal galaxies occupy a differ- 
ent locus in the colour-colour diagram. However, we note that 
the remaining ~ 30% of the 24 /im -detected X-ray-classified 
AGN/QSO also have the colours characterising normal galax- 
ies. Thus, we conclude that the application of a [5.8] - [8.0/im] 
colour cut can be useful to complement the X-ray AGN classi- 
fication, but cannot be both complete and reliable at the same 
time (cf. e.g. Barmby et al. 2006). 

In order to achieve a balance between completeness and re- 
liability for the active galaxy identifications within our sam- 
ple, we adopt the following criteria for our 24 yum -detected 
Ks < 21.5 mag sources. At redshifts z < 1.5, we base 
our AGN/QSO classification exclusively on the X-ray data. 
For sources at higher (z > 1.5) redshifts, we considered as 
AGN/QSO to all those objects classified as such in the X-rays, 
plus the two additional 24 yum sources which also have a colour 
[5.8] - [8.0/im] > 0.2 (AB mag; cf. figureE). In addition, we 
also classify as active galaxies eight additional objects of our 
sample which have been identified as AGN in the catalogue 



constructed by Alonso-Herrero et al. (120061 table 1 of this pa- 
per). Thus, using all these criteria, the total number of sources 
classified as AGN/QSO in our 24/im -detected Ks < 21.5 mag 
sample is 26 -H 2 + 8 = 36. Throughout this paper, when we 
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Table 1. The comoving number densities of massive galaxies with bolometric IR luminosity L/j? > 10" L© (LIRGs and ULIRGs) 
versus redshift. The redshift bins indicated with an asterisk are not complete for this luminosity cut, so the corresponding number 
densities should be taken as lower limits. The percentages given in the columns %{z) refer to the fraction of objects classified 
as LIRGs or ULIRGs among all the K^ < 21.5 mag galaxies in that redshift bin and with the same mass cut. The percentages 
%(z - 0) are obtained by dividing the number density of LIRGs and ULIRGs in that redshift/mass bin by the number density 
of galaxies with the same mass cut in the local Universe. These numbers indicate what fraction of the massive galaxies we see 
today have been IR-active in the past. Between brackets, we show the corrected values obtained when excluding the known 
AGN/QSO. These AGN/QSO are only those remaining in the 24 /im -detected Ks < 21.5 mag sub-sample considered here 
(i.e. with satisfactory HYPERZ SED fits and, thus, reliable mass estimates), which have been classified following the criteria 
explained in Section|3] Therefore, they constitute a lower Umit to the total number of mid-IR detected active galaxies. 
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Fig. 3. The comoving number densities of galaxies with bolometric IR luminosity L/s > lO" L© (filled symbols), with assembled 
stellar mass M > 1 x lO" M© (left panel) and M > 2.5 x lO" M© (right panel), versus redshift. Circles and upward-pointing 
triangles indicate directly derived values and lower limits, respectively. Both star-forming galaxies and AGN/QSO are considered 
in all cases. The empty circles correspond to the densities of the total K, < 21.5 mag galaxy population with the same mass cuts 
JC06a> . The asterisks at z = indicate the corresponding local values (Cole et al. I20(jll Bell et al. l2003> . 



refer to the AGN/QSO of our sample, we always consider all 
of these 36 objects, unless otherwise explicitly stated. 



4. The role of LIRGs and ULIRGs in stellar mass 
evolution 

We now consider our 24 //m -detected Ks < 2 1 .5 mag sample to 
study how the number density of massive luminous IR galax- 
ies evolved through cosmic time. Stellar masses have been ob- 
tained from the optical/near-IR SED of each galaxy, which in- 
cludes IRAC data at 3.6 and 4.5 |/m. The stellar masses are 
computed using the rest-frame A'^-band luminosities, which are 
interpolated/extrapolated from the closest observed passband. 
We note that the use of IRAC data guarantees that we directly 
map the near-IR light of galaxies up to high (z ~ 3) redshifts. 



We estimate our stellar masses to be accurate within a factor 
~ 2 (see IC06al for more details). 

Bolometric IR luminosities (L/s) for normal galaxies have 



been obtained using the Chary & Elbaz (1200 1> and Elbaz et 
al. (20021 formulae, as in C06b For the AGN/QSO present 
in our sample, we used the L(12yum) - Lir relation derived by 
Takeuchi et al. (2005 1, which makes no discrimination on the 
SED type. To interpolate/extrapolate the rest-frame 12/im flux 
of each AGN/QSO, we assumed that the SED of these objects 
followed a power-law fy cc V (a < 0). For each source, we 
derived the a value using its IRAC 3.6 and MIPS 24fj.ni fluxes. 
To assure that the luminosities obtained with the Takeuchi 
et al. ("2005^ relation were consistent with the Chary & Elbaz 
and Elbaz et al. formulae, we compared the obtained L/s values 
on our normal galaxies at z ~ 1, where the rest-frame I2fj.m 
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Fig. 4. Upper panel: the evolution of the stellar mass complete- 
ness with redshift, as imposed by the Ks = 21 .5 mag cut. These 
completeness limits are strict, as they have been computed us- 
ing a maximal mass-to-light ratio, corresponding to the tem- 
plate of a single stellar population formed at redshift ^ — > oo. 
Lower panel: the evolution of the bolometric IR luminosity 
completeness with redshift, as imposed by the completeness 
limits of the 24 jum survey. The luminosities in this plot have 
been computed using the Chary & Elbaz and Elbaz et al. for- 
mulae, and extrapolated at z < 0.4. 



flux is directly mapped by the observed 24 //m flux. We found 
that the luminosities obtained with the Takeuchi et al. relation 
are systematically lower by a factor ~ 1 .5 with respect to the 
luminosities obtained with the Chary & Elbaz and Elbaz et al. 
formulae. Thus, in order to maintain the same luminosity scale 
for all of our objects, we multiplied the AGN/QSO luminosities 
derived with the Takeuchi et al. relation by a factor 1.5. 

Figure |3l shows the redshift evolution of the comoving 
number densities of galaxies with bolometric IR luminos- 
ity L/R > lO" Lq (filled symbols) with assembled stellar 
mass M > 1 X 10" Mq and M > 2.5 x 10" Mq (left and 
right panels, respectively). By definition, the IR luminosity 
cut adopted for these galaxies implies that all of them are ei- 
ther LIRGs or ULIRGs. Both star-forming galaxies and ac- 
tive nuclei have been considered for all the number densi- 
ties shown in Figure |3 For star-forming galaxies, an IR lu- 
minosity L/R > 1O"(1O'^)L0 implies a star- formation rate 
SFR > 17(170)Moyr-i (Kennicutt 1998^. Circles in both pan- 
els of figure|3]correspond to directly derived values, i.e. number 
densities in the redshift bins in which we have completeness 



for all the L/« > lO" Lq galaxies (considering the 5v(24yum) - 
83/iJy limit). Upward-pointing triangles indicate lower limits 
on the densities in the redshift bins where our sample is not 
complete for the adopted luminosity cut. The evolution of the 
luminosity completeness imposed by the 5v.(24/im) = 83/iJy 
limit is shown in the lower panel of figure|3 For a comparison, 
the empty circles in both panels of Figure|3]show the densities 
of the total Ks < 21.5 mag sample, as they were computed by 
IC06a Analogous plots for only ULIRGs (Lm > 10'^ L©) are 
shown in Figure |5] 

Our computed number densities are ymflA:-corrected with re- 
spect to the Ks = 21.5 mag magnitude limit, to encompass 
the computation of number densities for the total Ks < 21.5 
mag galaxy population. No corresponding y„,fl (-corrections 
need to be applied to the 24 //m band, as we considered all the 
24 jum identifications in our sample with no flux cut. Instead, 
we use the 80% completeness limit of the 24 //m catalogue 
(Sv(24fj.m) = 83/vJy) as a reference to determine up to which 
redshift the 24 //m sample is complete for Ljr > 10" Lq and 
LiR > 10'^ Lq galaxies. We find that our 24yum sample is com- 
plete for LIRGs only up to redshift z ~ 1 and for ULIRGs up to 
redshift z ~ 2. Beyond those redshifts, we considered our com- 
puted number densities to be lower limits on the real values. 

We note that mass completeness is not a concern in this 
analysis. The upper panel of figurel^shows the evolution of the 
stellar mass completeness imposed by the Ks = 21.5 mag cut. 
For stellar masses M > 1 x 10" Mq, we are basically complete 
up to redshift z-3 (strictly M > 1.3 x 10" Mq at z = 3). 

The obtained number density values of Ljr > 10" L© and 
LiR > 10'^ Lq massive galaxies are listed in Tables ^ and |2 
respectively. In the last redshift bin z = 2-3, the depth of 
our 24 //m catalogue only allows us to detect ULIRGs, so the 
lower limits we obtain for the densities of Lm > 10" Lq and 
L/R > 10'^ Lq galaxies are the same at these redshifts. The 
results presented in both tables show that luminous IR galax- 
ies constitute a significant fraction of the assembled massive 
galaxies and this fraction increases with redshift. In particular, 
a minimum of ~ 65% of the most massive (M > 2.5 x 10" Mq) 
galaxies already present at redshifts z = 2-3 are ultra-luminous 
in the IR. 

Moreover, from comparison of figures |3] and [Sj we can 
clearly confirm how the role of IR galaxies among massive 
galaxies changes from being ULIRG-dominated at redshifts 
z > 1.5 to be LIRG-dominated at z < 1.5. In effect, we observe 
that the number density of massive ULIRGs drastically falls 
below z - 1-5. Within our sample, we find only one galaxy 
with LiR > IO'^Lq and M > 1 X 10" Mq at 1.0 < z < 1.5, and 
none below these redshifts. ULIRGs might be present at lower 
redshifts, but are indeed very rare (e.g. Flores et al. il999j . 

Within our sample, the ratio of AGN/QSO among the 
ULIRGs with stellar mass M > 2.5 x 10" Mq reaches its maxi- 
mum at z = 2 - 3. This fact, as well as the high S FR derived for 
the non AGN/QSO sources ( C06b i, highlight the importance 
of both star- formation and accretion activity in the evolution of 
the most massive galaxies already assembled at high redshifts. 
In contrast, the ratio of AGN/QSO among the Ljr > lO" Lq 
galaxies with stellar mass M > 1.0 x 10" Mq reaches its max- 
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Table 2. The comoving number densities of massive galaxies with bolometric IR luminosity L/r > 10'^ Lq (ULIRGs) versus 
redshift. See caption of Table^for references. 
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p,(xlO-^Mpc-3) 



%(z) 



% AGN(z) %(z = 0) 



p,(xlO-^Mpc-3) %(z) %AGN(z) %(z = 0) 



0.5-1.0 

1.0-1.5 0.58 ±0.58 (0.58 ±0.58) 0.8(0.8) 

1.5-2.0 15.0 ±2.7 (13.1 ±2.5) 35(30) 

2.0-3.0* 10.1 ± 1.5 ( 8.9 ± 1.4) 53 (47) 



0.8 (0.8) 0.58 ± 0.58 (0.58 ± 0.58) 0.8 (0.8) 

5 12(10) 7.3 ±1.9 (6.3 ±1.7) 47(41) 

6 8(7) 3.0 ±0.8 (2.3 ±0.7) 65(50) 
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Fig. 5. The comoving number densities of galaxies with bolometric IR luminosity L[r > lO'^L©. The references are the same 
as in Figure|3 In the left panel, the diamond at z = 2 shows the estimated density of ULIRGs obtained by Daddi et al. t.2005j . 
The cross-like symbol in the right panel indicates the number density of radio-detected sub-millimetre galaxies (Chapman et 
al.l2003l. 



imum at lower redshifts z - 0.5 - 1 .0, but is a minor fraction 
(6%) of all the .^5 < 21.5 mag galaxies with the same mass cut 
in the same redshift bin. 

With respect to the local population of galaxies with stel- 
lar mass M > 1.0 X 10" Mq (M > 2.5 x 10" Mq), we find 
that a minimum of 20%, 15%, 8% (40%, 36%, 15%) had IR 
luminosities Ljr > lO" Lq by redshifts z ~ 1.0, 1.5 and 2.0, 
respectively. 

Galaxies with Ljr > 10" L© contain ~ 24% of the total 
stellar mass density at redshifts z - 0.5 - 1 .0 and a minimum of 
~ 35 and 45% at z = 1. - 1.5 and 1.5 - 2.0, respectively. The 
total assembled stellar mass densities are 80-90%, 45-50% and 
25-30% of the local value at redshifts z ~ 0.75, 1.25 and 1.75, 
respectively (Caputi et al.[2005, C06a and references therein). 
This implies that L/r > 10" L© galaxies contain ~ 21%, and 
a minimum of 16 and 11% of the local stellar mass density at 
these respective redshifts. 

In the left panel of Figure |5] we compare our computed 
number densities of ULIRGs with the estimate obtained by 
Daddi et al. (2005 1 from the analysis of BzK-selected galax- 
ies (diamond-like symbol at z = 2), and we find a very good 
agreement. 

Finally, in the right panel of Figure |5l we show for ref- 
erence the number density of radio-detected sub-millimetre 



galaxies with flux 5y(850//m) > 5mJy (Chapman et al. 120031 
cross-like symbol in Figure |5}. The number density of these 
galaxies appears to be smaller by a factor ~ 4 (2-3) than the 
density of M > 2.5 x lO" Mq ULIRGs (L > 5 x IO'^Lq 
galaxies) at similar redshifts. There is evidence that the ma- 
jority of radio-detected sub-milUmetre galaxies have luminous 
24 /zm counterparts (Egami et al. l2004t and could be associated 
with the build up of the most massive galaxies (e.g. Stevens et 
al."2003'). However, the comparison of number densities shows 
that not all the most massive ULIRGs are associated with bright 
sub-millimetre galaxies (at least with radio-detected ones). 
What kind of 24/im -selected ULIRGs are counterparts to sub- 
millimetre galaxies is still not completely clear (cf. Lutz et 
al. I2005> . and this is being investigated, for instance, as a part 
of the SCUBA Half Degree Extragalactic Survey (SHADES) 
project (Mortier et al. 120051 



5. Linking the ULIRG and the ERG phases 

The colours characterizing ERGs (/ - Ks > 4, Vega mag) are 
usually associated with the presence of old stellar populations 
and/or starbursts heavily enshrouded by dust. 24 //m data can 
be used to constrain the importance of these two components 
in the ERG population. 
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Fig. 6. The (/ - Ks) colours versus redshift for the Ks < 21.5 mag galaxies in different stellar mass bins. The /-band data 
corresponds to ACS/HST data. The filled symbols indicate 24yum -detected ULIRGs. All the triangles correspond to lower limits 
on the (/ - Ks) colours for sources below the ~ 2cr detection limit in the / band. In each panel, the dotted vertical line shows the 
strict completeness redshift limit of the corresponding stellar mass bin. 



Figure |6l shows the (; - Ks) colours of the Ks < 21.5 
mag galaxies of different stellar mass, with the 24/im -detected 
ULIRGs indicated by filled symbols. Although formally the 
Chary & Elbaz (1200 It and Elbaz et al. (EUU^ relations do not 
allow us to compute bolometric IR luminosities beyond z ~ 3, 
the limits of the 24 fim survey show that only ULIRGs can be 
detected beyond z ~ 2. Thus, every 24 yum -detected source 
at z > 3 is indicated as a ULIRG in Figure |6l Following our 
adopted classification criteria (cf. Section|3}, about one third of 
our 24 yum -detected galaxies at z > 3 are AGN/QSO. However, 
we note that our classification might be particularly incomplete 
at these high redshifts and the fraction of active galaxies among 
the 24 yum -detected galaxies might be higher. Unfortunately, 
the present data do not allow us to better constrain the AGN 
fraction. 

We observe that ULIRGs are associated with ^5 galaxies 
with a variety of stellar masses and (/ - Ks) colours at redshifts 
z> 1.5, with ~ 60% of them having the characteristic colours 
ofERGs. 

Conversely, at 1.5 < z < 2, ~ 29% of the Ks < 21.5 mag 
ERGs are ULIRGs. At 2 < z < 3, this percentage increases to a 
minimum of 49%. If we consider only those ERGs with stellar 
mass M > 1.3x10" Mq, these fractions are 40% at 1 .5 < z < 2 
and > 52% at 2 < z < 3. In particular, from Figure|6lwe see that 
some of the reddest ERGs are not ultra-luminous in the IR (and 
some are not detected at all in our 24 yum sample). Put another 



way, no correlation seems to exist between the (/ - Ks) colours 
and the IR luminosity. This result appears to be independent of 
the assembled stellar mass. We note that the identification of 
ERGs with ULIRGs is strictly complete up to redshift z ~ 2. 
Above that redshift, our 24yum sample starts to lose complete- 
ness even for ULIRGs, so there exists the possibility that some 
of the non-24 //m detected ERGs are also ultra-luminous in the 
IR. The presence of old stellar populations might in part be re- 
sponsible for the (/ - Ks) colour excess of those ERGs that are 
not ULIRGs. However, it is known that many of the reddest 
colours observed in ERGs can only be explained by the super- 
position of an evolved stellar component and non-negligible 
amounts of dust (Caputi et al. l2004l Papovich et al. l2006> . 



Daddi et al. J2005> find the average duty cycle for vigor- 
ous starburst with stellar mass M > 10" Mq to be at least 50% 
at redshifts 1.4 < z < 2.5. In this context, the duty cycle is 
the fraction of time spent in the starburst phase over the to- 
tal elapsed time in the redshift interval 1.4 < z < 2.5. The 
percentages of ULIRGs among massive ERGs we find here are 
consistent with the duty cycle estimated by Daddi et al.. Recent 
studies determined that most of the brightest 24 jum galaxies at 
z > 2 akeady have stellar masses M > 1 x 10" Mq and that 
star-formation activity should mainly have proceeded by mul- 
tiple discrete starbursts ('C06b' Papovich et al. I2006> . On the 
other hand, the majority of the massive galaxies at z > 1 .5 dis- 
play extremely red colours. All these facts converge to a see- 
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nario in which most of the ERGs -and not only the 24/im- 
detected ones- could suffer a bright 'mid-IR episode' (i.e. a 
ULIRG phase) up to several times in their lives. After a pe- 
riod of star-formation activity, a galaxy becomes fainter in the 
mid-IR, but it still would be an ERG. 

We compared the SEDs of our ERGs identified as ULIRGs, 
with those belonging to the other ERGs with similar stellar 
masses and redshifts. We found that the average properties 
(star formation history, age) of these SEDs are the same in the 
two groups of objects. Even the colour excesses are consistent 
within the error bars. The colour excesses are obtained from the 
convolution of the template SED with the Calzetti et al. ( 2000j 
reddening law. For example, if we consider those ERGs with 
stellar mass M > 1.3 x lO" Mq at redshifts 2 < z < 3 that 
are ULIRGs and that are not, we find that the medians of the 
colour excess of the best-fit SEDs are E{B - V) = 0.30 + 0.07 
and E(B - V) = 0.22 + 0.05, respectively. For the ERGs with 
the same mass cut but lying at redshifts 1 < z < 2, the respec- 
tive medians of the colour excess are E{B - V) = 0.17 + 0.05 
and E{B - V) - 0.10 + 0.07. These results are basically 
the same if we divide the ERGs into those which have flux 
5v(24//m) > 83 //Jy and 5y(24/^m) < 83 //Jy (including, of 
course, 24 yum non-detections). 

Yan et al. (2004i suggested that the 24 //m -detected ERGs 
could be in the process of transforming into non-24 /im- 
detected early-type ERGs at z ~ 1. Our results are consis- 
tent with such a transformation hypothesis for the ERGs. More 
specifically: 1) even if evolved stellar populations were already 
present at 2 < z < 3, the E{B - V) colour excesses character- 
ising massive ERGs at these redshifts would be inconsistent 
with these objects being passive systems. It is possible that 
most of these objects are actually IR luminous, and some of 
them are non-detected within the limits of our survey. However, 
if star formation indeed proceeds by discrete episodes, it will 
be more likely that a substantial fraction of the non-24 //m- 
detected ERGs are the detected ones in a post-starburst phase. 
2) at 1 < z < 2, the colour excesses are smaller on average, 
indicating that old stellar populations have a more important 
role in the colours of ERGs. However, some of the ERGs at 
these redshifts could still correspond to post-starbursts, as it 
has been demonstrated by spectroscopic studies (e.g. Doherty 
etal. l2005l . 



6. Summary and discussion 

In this work we quantified the importance of luminous IR 
galaxies (LIRGs and ULIRGs) in the evolution of the massive 
A'j-band galaxy population. This allowed us to probe the role 
of star-forming galaxies (and quasar activity) in the history of 
stellar mass assembly. 

Our results show that a substantial fraction of the massive 
galaxies (M > 1 x 10" Mq) assembled at diff'erent redshifts are 
luminous in the IR. This fraction appears to be higher at high 
redshifts (z ~ 2-3)thanatlowerredshifts. This resultis clearly 
indicating that IR activity in massive galaxies -either due to 
star formation or AGN-driven- has been much more important 
in the past. 



We determined that a minimum of 53 and 65% of the most 
massive (M > IxlO'' MoandM > 2.5x10" Mq, respectively) 
galaxies present at z ~ 2 - 3 were ULIRGs at those redshifts. 

Although the formation redshift of the first most massive 
systems cannot be determined from the depth of current sur- 
veys, our results show that the elapsed time between z = 3 and 
z = 2 is of major importance in the construction of the bulk 
of the stellar mass in a significant fraction of massive galax- 
ies. The number density of all the galaxies with stellar mass 
M > 1 X 10" Mq (M > 2.5 x 10" Mq) rises by more than a 
factor two (three) between z ~ 3 and z ~ 1.5 (cf. fig.|3and|5}- 
This is entirely consistent with a major period of IR activity at 
z~2-3. 

The similarities between the optical/near-IR SEDs of mas- 
sive ERGs that are ultra-luminous in the IR and those that are 
not suggests that ULIRGs do not represent a particular type of 
near-IR selected galaxies. They rather could constitute a phase 
(the phase of high star-formation/quasar activity) during the 
life of a massive galaxy, and in particular during the period 
in which it is detected as an ERG. Bright IR galaxies are only 
able to trace a fraction of the massive galaxies present in the 
Universe at a given time, but this fraction becomes very impor- 
tant ( > 50%) at redshifts z > 2. 
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